Anisotropic nonlinear elasticity in a spherical bead pack: influence of the fabric anisotropy 
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Stress-strain measurements and ultrasound propagation experiments in glass bead packs have been simulta- 
neously conducted to characterize the stress-induced anisotropy under uniaxial loading. These measurements, 
realized respectively with finite and incremental deformations of the granular assembly, are analyzed within the 
framework of the effective medium theory based on the Hertz-Mindlin contact theory. Our work shows that 
both compressional and shear wave velocities and consequently the incremental elastic moduli agree fairly well 
with the effective medium model by Johnson et al. [J. Appl. Mech. 65, 380 (1998)], but the anisotropic stress 
ratio resulting from finite deformation does not at all. As indicated by numerical simulations, the discrepancy 
may arise from the fact that the model doesn't properly allow the grains to relax from the affine motion ap- 
proximation. Here we find that the interaction nature at the grain contact could also play a crucial role for the 
relevant prediction by the model; indeed, such discrepancy can be significantly reduced if the frictional resis- 
tance between grains is removed. Another main experimental finding is the influence of the inherent anisotropy 
of granular packs, realized by different protocols of the sample preparation. Our results reveal that compres- 
sional waves are more sensitive to the stress-induced anisotropy, whereas the shear waves are more sensitive to 
the fabric anisotropy, not being accounted in analytical effective medium models. 
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PACS numbers: 45.70. -n, 43.35.+d, 81.05.Rm 
I. INTRODUCTION 

The study of nonlinear elasticity in granular materials is not 
only of fundamental interest but also of practical importance 
in many fields, including soil mechanics and geophysics. Un- 
like consolidated porous materials such as sedimentary rocks, 
soils and sintered bead packs, dry granular materials acquire 
solely its elasticity as a result of the applied stress, forming 
a very inhomogeneous force network yl [2, [H 0, [H 0] . The 
mechanical properties of such materials depend strongly on 
both the contact interactions between grains and the geomet- 
ric arrangement (i.e. packing structure). The stress depen- 
dence of elasticity which originates with these contact forces 
at the grain level leads to a complex behaviour of the granular 
medium: nonlinear behaviour, loading-path dependence and 
stress-induced anisotropy 0]. The change of sound velocity 
with the applied stress (i.e. acousto-elastic effect)is the char- 
acteristic signature of this nonlinear elasticity HI 0, [H [H] . The 
hysteresis often found in the stress-strain experiments of gran- 
ular media is also associated with the path-dependant nature 
of the contact forces Jfjl [HI EH] ■ Another issue complicating 
the analysis of the mechanical properties of granular materials 
arises from the geometric or fabric anisotropy, which is basi- 
cally related to the mode of grain deposition and the packing 
structure in a model system like the sphere pack ifiol MJl. This 
fabric anisotropy is important for the mechanical behaviour of 
granular media such as stress transmission lfl3ll . mechanical 
stability iM and liquefaction resistance lfl5ll . 

Photoelastic visualizations lfl6ll and numerical simulations 
in 2D systems lfl7l[T8l[l9il20ll demonstrate that the anisotropy 
induced by an external load may have two distinct effects: one 



leads to an important change of the fabric anisotropy in the 
contact network, and the other develops an anisotropic force 
chain network. Understanding the nonlinear elastic responses 
and the associated stress-induced anisotropy represents a fun- 
damental issue for granular mechanics. The problem may also 
be important for understanding the jamming phenomena and 
shear-induced yielding in a more general class of geomateri- 
als formed from granular media IU4l[l5l[l6ll . Development of 
non destructive methods of investigation is therefore desirable 
for monitoring the evolution of both geometric and mechani- 
cal anisotropies in real 3D granular materials, and allows us to 
gain a more comprehensive insight into their elastic properties 
along different stress paths. 

Sound waves offer a sensitive and non invasive probe of 
both the structure and the mechanical properties of hetero- 
geneous materials ll2ll I22T l23i l24ll . In a granular medium, 
sound propagation is controlled to a large extent by the prop- 
erties of the contact force networks forming the solid frame of 
the material 0, H HHH H3] . The effective medium theory 
(EMT) has been commonly used to describe both the nonlin- 
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and sound propagation in the 
I3dl3lll . These effective medium 
models are generally based on the Hertz-Mindlin theory at the 
grain contact and make use of the two main assumptions to 
obtain a mean field description of the granular elasticity: i) 
affine approximation or kinematic hypothesis, in which the 
motion of each grain follows the applied strain, and ii) sta- 
tistically isotropic distribution of contacts around each grain. 
However, some acoustic measurements in random packs of 
glass beads under isotropic loading j4l [32Tl show that the ratio 
of bulk modulus K to shear one G is significantly larger than 
the value predicted by the EMT with bonded elastic spheres 
(no sliding). Numerical simulations reveal that the discrep- 
ancies can stem from the failure of the effective medium ap- 
proximation: the bulk modulus is well described by the EMT 
but the shear modulus is not, principally because the EMT 
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does not correctly allow the grains to relax collectively from 
the affine motion assumed by the theory 0,[33t]- In a recent 
work, the multiple scattering of shear acoustic waves through 
the stressed glass bead packs shows that sound propagation at 
small-amplitude (or incremental deformation) does not cause 
any significant rearrangement of the contact force network 
12a]. This observation leads to a fundamental question as to 
whether overestimating the shear modulus by the EMT is re- 
lated with the failure of the effective medium approximation 
for sound propagation or with the inadequate treatment of mi- 
cro mechanics at the grain contact level 171, 12211 . 

In this paper, we present new results obtained from a granu- 
lar model system, i.e. a glass bead pack in an cedometric test. 
We conducted simultaneously the stress-strain measurement 
and the ultrasound propagation experiment to characterize the 
stress-induced anisotropy by uniaxial loading. Both measure- 
ments of stress fields and elastic moduli via sound velocities 
parallel and perpendicular to the applied load are analyzed 
within the framework of the EMT developed by Johnson et 
al. specifically for stress induced anisotropy in the cedometric 
configuration l29ll . The aim of our work is twofold. Firstly, 
we examine the applicability of the EMT to small- and large- 
deformation mechanical tests, corresponding respectively to 
sound propagation and stress-strain measurements. The cru- 
cial role of the interaction nature at the grain contact for the 
appropriate prediction will be discussed. Then, we investigate 
the respective responses of compressional and shear waves to 
the elastic anisotropy in granular media. A particular atten- 
tion is paid to the influence of the fabric anisotropy which is 
realized by the different sample preparation. Such inherent 
fabric anisotropy is expected to evolve little under cedometric 
loading compared to pure shear experiments. 



II. EXPERIMENTS 

In order to characterize the anisotropic elasticity of granular 
materials, we developed an apparatus coupling the mechanical 
test with the ultrasonic measurement. The schematic diagram 
of the experimental arrangement is illustrated in Fig. Q] Spher- 
ical glass beads of diameter between 300 jim and 400 /mi are 
carefully filled in a duralumin cell up to 70 mm in height and 
closed by the fitted pistons at the top and bottom surfaces. 
A controlled vertical load is applied to the granular sample 
across the top piston. The forces re-directed in the horizontal 
direction are measured by the stiff force sensors placed at the 
middle of the cell to monitor the stress anisotropy evolution. 
This special compaction cell is presumed to be rigid enough 
that the sample experiences a macroscopic uniaxial strain with 
negligible lateral motion (cedometric test). The solid volume 
fraction <p of the granular assembly is determined by the axial 
displacement as a function of the applied load. To ensure the 
measurements of sound velocities in both the vertical (z-axis) 
and horizontal (x-axis) directions, a pair of compressional or 
shear piezoelectric transducers of diameter 30 mm in contact 
with the glass beads are placed at the top and bottom pistons, 
while the other ones are mounted at the lateral walls separated 
by a distance of 40 mm (see the inset of Fig. [TJ. 
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FIG. 1: Sketch of the experimental set-up. Inset: Top- view of the 
cedometric cell. 



A. Granular fabric and packing density 

To study the influence of the granular fabric on the induced 
elastic anisotropy, the bead packs are realized by two different 
preparation protocols. The first one, known as rain deposition 
or air pluviation, consists in pouring the glass beads into the 
cell through two grids of mesh size 2 mm spaced by 2 cm and 
reaches a dense pack of <p ~ 0.64. The second, termed as 
de-compaction protocol is to remove gently a horizontal grid 
through the bead pack from the bottom to the top after fill- 
ing the beads in the cell; the rearrangements produced by this 
efficient shearing throughout the sample allow us to obtain a 
loose pack of <p * 0.60. It is shown by numerical simulations 
|G3 H that the rain deposition protocol creates an anisotropic 
distribution of the contact angle with two preferred directions 
orientated roughly at 30° around the vertical (gravity) direc- 
tion. However, the loose pack prepared by the de-compaction 
protocol is expected to p roduce a fairly isotropic distribution 
of the contact angle 13511 . 

Before any measurement, a preloading up to cr zz = 400 kPa 
is applied to the sample in order to minimize the hysteretic 
effects related to the grain rearrangements and ensure a re- 
producible initial state. Then the stress-strain and ultrasonic 
velocity measurements are performed as a function of the ap- 
plied stress cr zz ranging from 70 to 900 kPa. Figure|2]displays 
the evolution of the solid volume fraction (p versus cr zz in the 
dense and loose packs, obtained respectively by the two prepa- 
ration protocols. These results show a good reproducibility of 
our sample preparation in terms of the initial packing den- 
sity measured at cr zz = 70 kPa for 30 repeated measurements: 
(pdense = 0.642 ± 0.002 and </>i oose = 0.605 + 0.002. The fact 
that the solid volume fraction in the dense packing sample is 
a little larger than the density of the random close packing 
(RCP) is due to the slight dispersion of the bead size. Fur- 
thermore, Fig. |2] shows that the packing density of our sam- 
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FIG. 2: Evolution of solid volume fraction of the loose and dense 
packs. Error bars illustrate the data dispersion over 30 measurements. 
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FIG. 3: Evolution of the vertical-to-horizontal stress ratio cr^/cr^ 
versus the applied stress cr- in the dense and loose packs. 



pies increases slightly, up to a resultant vertical deformation 
of sj= e 3 ) = 3 x 1(T 3 at the applied (r a (= cr 3 ) = 900 kPa. 
Both the dense and loose bead packs keep clearly a signature 
of their preparation protocols. These results imply that for the 
applied stress in this work, much less than the value to pro- 
duce the grain fracture of about 20 MPa for the glass beads 
l36ll . there are few important rearrangements of grains in the 
cedometric test as compared to those in a pure shear exper- 
iment Hill and the inherent fabric should not change much 
from the initial state 



B. Anisotropic stress field 

Let us now investigate the stress-induced anisotropy under 
uniaxial loading in the cedometric configuration. By moni- 
toring the resultant horizontal stress cr xx when increasing the 
vertical stress <t zz , we display in Fig. [3]the stress ratio cr zz /cr xx 
measured as a function of the applied stress <r zz for the two 
granular samples, respectively. It is observed that both the 
bead packs evolve from their initial stress fields at low ap- 
plied cr zz = 70 kPa, roughly isotropic, to anisotropic states 
at high cr zz . As mentioned in the previous experiment l29ll . 
the isotropic distribution of stress observed at low cr zz arises 
probably from a kind of tight wedging of grains produced dur- 
ing the preloading and unloading cycle. Such an interlock 
of grains may contribute to a residual stress isotropically dis- 
tributed in the bead pack even when the top piston is removed 
and no stress is applied. 

As the applied stress increases, the ratio cr zz /cr xx is seen to 
increase to an asymptotic value at high cr zz . For the granu- 
lar packs prepared by the two distinct protocols, the different 
asymptotic values of about 1 1 % reveal again a memory effect 
of the initial state of the sample as seen in the above density 
measurement. In the conventional cedometric test, the inverse 
of the asymptotic value cr zz /o- xx is known as Jacky coefficient 
Kq of the earth pressure at rest; it is empirically related to the 



internal friction angle 6q of the medium by ,Ko = 1 - sin 0q Ir37ll . 
These results thus suggest a possible correlation between the 
internal friction angle and the fabric of the medium. To ex- 
amine this issue, we have measured the avalanche angle 6 at 
the relative humidity RH = 40% for both the dense and loose 
packs, and obtained 9\ oose = 30.4+1.0° and = 33.2+0.3° 
for six avalanche experiments. The difference between the 
two measured avalanche angles is about 9%, consistent with 
the difference of 11% between the two internal friction an- 
gles deduced from the asymptotic values of cr zz /cr xx in Fig. 
[3] Although the avalanche angle is not the direct measure of 
the internal friction angle, these stress measurements indicates 
again the significant effect of the preparation protocol on the 
granular fabric 03. 



C. Anisotropic elastic moduli 

Different transports of sound waves in granular media have 
been detailed in the previous work lf39ll . When the wavelength 
is much larger than the bead size, the coherent waves propa- 
gate as through an effectively homogeneous medium. Mea- 
suring the elastic wave velocity V allows one to access to the 
elastic modulus M of the granular medium by M = pV 2 with p 
the material density m^ESl- In our apparatus, the compres- 
sional P-waves and shear S-waves are excited and detected, 
respectively, by longitudinal or transversal piezoelectric trans- 
ducers of diameter 30 mm. Figure|4]displays a typical pulsed 
ultrasound transmission along the vertical direction. The ar- 
rivals of P- and S-wave pulses are well separated, which al- 
lows us to clearly identify the different modes and measure 
adequately the wave velocities by the time-of-flight. 

For the granular materials under uniaxial loading and the 
geometry of the cell considered here, the elastic tensor C/j 
relating the incremental stress dcri and strain dsj, dcrj = 
Cjjdsj, would have a symmetry belonging to the orthotropic 
class. In such a case, there are nine independent constants as 
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FIG. 4: Typical ultrasound transmission along the vertical direction 
(<x. ; =500 kPa), excited and dectected by shear transducers. 
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By means of P- and S-waves of different polarization propa- 
gating along the horizontal and vertical directions (i.e. princi- 
pal axes), respectively, we can infer the diagonal elastic com- 
ponents (/ = J) from the velocity measurement C// = pVfj. 
For determining the off-diagonal components, velocity mea- 
surements of waves propagating along a direction inclined to 
the principal axes should be necessary. We display in Fig. 
|5h and |5]3 the wave velocities of several elastic waves mea- 
sured a function of the applied stress o~ zz in our loose pack 
and dense one, respectively. Here V\ \ and V33 correspond to 
the velocities of P-waves propagating along the horizontal and 
the vertical (x- and z-axis) directions. V44 denotes to the ve- 
locity of the S-wave propagating vertically, whereas V55 and 
V(,(, correspond to those of the S-waves propagating horizon- 
tally and being polarized along the vertical and perpendicular 
(y-axis) directions. Each data presented in Fig. |5]results from 
an average of six experimental runs and the data dispersion or 
error-bar is lower than 2% in the case of the loose pack and 3% 
in the dense sample, illustrating thus a good reproducibility of 
these acoustic measurements. 

As shown in Fig. [5^ for the loose pack, there is a clear dif- 
ference of P-wave velocity between V\ \ and V33 (> V\\) going 
up to 15% when the stress cr zz is increased to 900 kPa. For the 
S-wave, the velocity difference V44 > V55 > is less im- 
portant than that of P-wave: indeed, at cr zz - 900 kPa V44 is 
about 5% and 8% greater than V55 and V66, respectively. In 
contrast to the previous works under load superior to a few 
MPa |0, HH [32[| . our acoustic measurements are realized at 
much lower stress, ranging from 70 kPa to 900 kPa. In such 
range of loading, it is expected that both the solid volume frac- 
tion <p (Fig. [2]i and the coordination number Z would vary little 



|01 ■ As the de-compaction protocol tends to create an isotropic 
fabric in the loose granular sample, our measurements indi- 
cate that the elastic anisotropy observed here shall stem from 
the induced stress anisotropy shown in Fig. [3] 

For the dense pack prepared by rain deposition, Fig. 
displays again a significant difference between the P-wave ve- 
locities Vu and V33. When the anisotropy of the stress field is 
developed, such elastic anisotropy rises to 18% at cr :z = 900 
kPa. Unlike the P-wave velocity, the S-wave velocity appears 
to be sensitive to the fabric anisotropy produced by rain de- 
position in the dense granular pack. Contrary to the mea- 
surement in the loose packing, we observe a significant dif- 
ference between V44 = 306 + 2 m/s, V55 = 342 + 7 m/s and 
V66 = 3 3 3 + 6 m/s at the lowest stress cr zz = 72 kPa where the 
stress field is nearly isotropic (see Fig. [3j. Moreover, when 
cr zz is increased the S-wave velocities evolve differently, but 
recover a similar behaviour at high load: a difference of 8% 
is found again between V44 and Vgg at cr zz = 900 kPa. These 
results imply a particular sensitivity of S-wave to the fabric 
anisotropy which affects the elastic anisotropy together with 
the stress anisotropy. 



III. COMPARISON WITH THE EFFECTIVE MEDIUM 
THEORY 

The theory of elasticity a granular pack is primarily based 
on the Hertz-Mindlin model of contact between grains |42il . 
The macroscopic stress-strain relations are commonly derived 
using the effective medium description, where it is assumed 
that the motion of grains is affine with the applied macro- 
scopic strain at least on average, and the distribution of con- 
tacts is statistically isotropic and homogeneous. Because the 
presence of the tangential forces at the contacts gives rise to 
load-displacement relations which are not only nonlinear but 
also inelastic, the mechanical response of the medium, namely 
stress-strain relations, must be expected to depend on the en- 
tire past history of loading. However, it is shown that the 
incremental response of the medium, i.e. the second-order 
effective elastic constants and consequently sound velocities 
are path-independent of loading (0]. Within the framework 
of the effective medium approach, Walton J28ll analyzed the 
mechanical responses of a random pack of identical elastic 
spheres under isotropic strain and purely uniaxial compres- 
sion, respectively. Furthermore the incremental elastic mod- 
uli were derived for these specific initial deformed states. For 
simplicity the spheres were assumed to be either infinitely 
rough or perfectly smooth in the calculations; nevertheless 
these results provide a physical insight of the effect of fric- 
tion in real grain contacts. 

As noted above, our measurements were made in the com- 
paction cell with rigid walls in which the granular sample is 
subjected neither to purely isotropic nor uniaxial compres- 
sion. To compare with the experimental data obtained in such 
widely used oedometric tests, Johnson et al. proposed an an- 
alytical model based on the effective medium theory I29I1 . 
This model combines isotropic and uniaxial strains to describe 
the stress-induced elastic anisotropy in transversely isotropic 
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FIG. 5: Wave velocity versus applied stress cr- in the loose pack (a) 
and the dense one (b). Experimental data (♦ : V 33 , • : V n , ■ : V 44 , 
▼ : Vj5 and A : V(,(,) are compared to the EMT predictions (lines) . 



for purely uniaxial compression in the opposite limit s — > 0, 
respectively ll28tl . 

Fig. [3] presents the ratio cr a lcr xx computed for the packs 
of both bonded ("rough") and frictionless ("smooth") spheres, 
as a function of the applied stress cr zz . The latter case is 
obtained by canceling the tangential contact stiffness in the 
model, i.e. C, — > 0. As described in eqs. 21 and 22 of 



029I1 . the stress tensor cr^ is a function of several parame- 
ters, (Tij = (Tij(s, S3, Z, <p, jj. g ,v g ) where Z is the coordination 
number, fj. g (= 24MPa) is the shear modulus and v g (= 0.2) is 
the Poisson ratio of the glass bead. For a given C, the ratio 
cr zz /(T xx versus the applied stress <t zz is only parameterized by 
the glass bead property fi g and v g (Fig. [3]). The axial strain 
£3 is measured from the displacement of the top piston ver- 
sus <t zz and the residual isotropic strain s is deduced from the 
isotropic residual stress 70 kPa measured at £3 = 0, yielding 
e= 1.5 X 10~ 4 . 

We observe from Fig. [3] that the theoretical curve calcu- 
lated with the frictionless spheres (C, — > 0) is faithful to the 
qualitative trend observed in our experiments and is overall in 
agreement with the measured data; at the high load the stress 
ratio approaches to the limit value of uniaxial loading, i.e. 
cr zz/ cr xx = 4 J28ll29ll . In contrast, Fig. [3] and the inset show 
that the no-slip assumption with bonded spheres gives rise to 
a huge overestimation of this stress ratio, indicating the cru- 
cial role of the contact interaction law in the effective medium 
model for describing adequately the evolution of the stress 
field anisotropy. We may understand the discrepancy found 
with the bonded spheres by the following picture. At large de- 
formation applied here of the order of 10~ 3 , the grains need to 
rearrange in order to relax the tangential stress accumulated at 
the contacts. The no-slip assumption with infinitely frictional 
spheres forbidden such a process and hence fail, whereas the 
effective medium analysis with frictionless spheres allow the 
stress relaxation and provide thus an adequate description of 
the experiments. 

Let us now examine the applicability of the effective 
medium approach to the wave velocity measurements. Given 
a stress state cr and corresponding macroscopic deforma- 
tion £, the elastic moduli can be derived from the incremen- 
tal stress of the medium subject to an incremental strain by 
do-jj = Cijktdski. Following the work of l29ll . yields 

3<pZ 



Qjki 



+B {8i k njni + Sun jn k + djinm + (J^n/n/))) ( 1 ) 



where B 



granular materials: Sij = sdij + £36136 j3 with 6ij = 1 if i = j, 
otherwise 5jj = 0. Here the internal strain £ stems from the 
residual stress distribution, namely an initial isotropic stress 
state produced by the preloading-unloading process as men- 
tioned above, whereas the axial strain £3 results from the ap- 
plied stress cr zz . As V55 differs from V44 only by 5% (Fig. 
0, the approximation of the transverse isotropy used in this 
model would be adequate for interpreting our measurements. 
The predictions of such a model reduce to the results for 
purely isotropic compression in the limit £3 — > and to those 
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- -xv - e nR, and <5, , denotes the 
Kronecker symbol. The brackets () represent an average over 
all unit vector n uniformly distributed since the distribution of 
contact angle (i.e. granular fabric) is assumed to be isotropic. 
To determine the wave velocities from the elastic moduli in 
eq. [T] we perform the calculations by two steps. Firstly, 
we compute the stress state of a frictionless spheres packing 
(C, = 0) created by the large axial deformation, which pro- 
vides an adequate description of the experimental stress field 
shown in Fig. [3] Then, we consider the appropriate value of C, 
{+ 0) according the Hert-Mindlin theory to calculate the elas- 
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tic moduli by eq. Q] Note that this procedure of computation 
is also used in many numerical simulations where the friction 
between the particles is turned off during the packing prepa- 
ration and is turned on to measure the incremental response 
01]. Fig. [5^ and[5t> illustrate the calculated wave velocities 
in both the loose packing and the dense one. The agreement 
between the theoretical predictions and experimental results 
remains fairly good (better than 10%), for the present effec- 
tive medium model employing only one adjustable parameter, 
i.e. Z = 6 in the dense pack and Z = 5.5 in the loose pack. 



IV. DISCUSSIONS 

A. Stress dependence of P-wave velocities 

Based on the Hertz-Mindlin theory of contact and the as- 
sumption of isotropic fabric in the sphere pack, the effective 
medium theory (EMT) predicts a power-law dependence of 
the elastic modulus M ~ P 1 ^ (or V ~ P 1 ' 6 ) on isotropic 
compression P 01 0, H3, HH]. The same power law scaling 
is found for the pure uniaxial loading (i.e. ey = esSqSjs) 

1/3 1/3 
as M ~ cr-l or M ~ cr xx since the stress ratio cr zz /cr xx is 

constant |[28tl . Numerous acoustic velocity measurements in 
sands and glass bead packs showed however that the exponent 
in the power-law scaling may not be constant, var ying from 
1 1 A at low pressure to 1 /6 at high pressure Illlsllil "Hi. Sev- 
eral mechanisms have been proposed to explain this discrep- 
ancy, including the contact recruitment by buckling of particle 
chains or the conical contact due to irregular surfaces 12], the 
effect of soft shell on coated spheres 14511 and the fluctuation 
of the stress field l46ll . However, for the preloaded glass bead 
packs and the range of the applied stress considered here, the 
contact recruitment is expected to have a negligible effect, as 
that suggested by the above density measurement (Fig. |2|. 

We investigate here the influence of the stress anisotropy on 
the scaling behaviour of the P-wave velocity. To do this, we 
displays in Fig. |6]the rescaled P-wave velocities V23 / all 6 and 
Vu/crll 6 as a function of the applied axial stress cr zz , obtained 
in the loose and dense packs, respectively. For a zz > 300 kPa, 
V33 is seen to follows adequately the scaling of a\ z 6 , as ex- 
pected within the framework of the Hertz-Mindlin theory. By 
contrast, V\ 1 deviates drastically from the all 6 scaling, which 
is due to the evolution of stress-field anisotropy illustrated in 
Fig. [3] Indeed, if rescaling Vu with the stress <t xx parallel 
to the propagation direction, we recover a scaling behaviour 
V11 ~ o"x( 6 (insets to Figs. [6^ and|6j5), similar to those for 
V33 ~ all 6 in both the loose and the dense packs having the 
different granular fabric. 



B. Stress and fabric dependences of S-wave velocities 

Compared to the P-wave velocity, the stress dependence of 
the S-wave velocity is more complicated and strongly influ- 
enced by the fabric of the pack. As mentioned above (in sec. 
Hit . S-waves are not only sensitive to the stress field anisotropy 
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FIG. 6: Stress dependence of P-wave velocities propagating verti- 
cally (♦ : V33) and horizontally (• : Vi\) in the loose pack (a) and 
the dense one (b). The dashed lines correspond to the predictions by 
Hertz-Mindlin theory of contact. 



but also to the fabric anisotropy. Let us first examine the scal- 
ing behaviour of the shear velocities versus a zz in the loose 
bead pack with isotropic fabric. As shown in Fig. [7^, both V44 
and V55 scale roughly as a power law i.e. ~ o%/ 6 for cr zz > 300 
kPa, which agree fairly well with the stress-dependence pre- 
dicted by the Hertz-Mindlin theory of contact. Compared to 
the P-wave velocity Vn, the scaling behaviours of the shear 
velocity V55 is somewhat surprising, however, it may be un- 
derstood by a heuristic picture which accounts for both the 
propagation direction and the polarization. Indeed, unlike the 
P-wave propagating in the horizontal direction (Vn), the S- 
wave travelling along the x-axis is polarized vertically along 
the z-axis (V55) and it would also be affected by the stress com- 
ponent cr zz . This scenario is consistent with the preceding ob- 
servation showing that the S-wave velocity is less sensitive to 
the stress-field anisotropy than the P-wave velocity is (shown 
in Fig. [5]). However, for the S-wave travelling along the x- 
axis but polarized horizontally along the y-axis (V66), Fig- Eh 



shows a scaling behaviour of V66 deviated from all 6 but being 
consistent with the power law ~ cr'{ 6 (inset to Fig. [7^), a 
situation reminiscent to that of the P-wave (Vn). 
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FIG. 7: Stress dependence of S-wave velocities propagating ver- 
tically (■ : V 44 ) and horizontally (t : V 55 and ▲ : V<&) in the loose 
pack (a) and the dense one (b). The dashed lines correspond to the 
predictions by the Hertz-Mindlin contact theory. 



We now turn on the scaling behaviour of the S-wave ve- 
locities in the dense pack presenting the fabric anisotropy ob- 
tained by rain deposition. Contrary to the behaviours observed 
in the loose packing with isotropic fabric, there is no clear 
stress dependence of the shear velocities V44, V55, and in 
the range of the applied stress (Fig. [7J5 and the inset). This 
observation reveals the extreme sensitivity of the shear wave 
velocity to the fabric of the granular pack. However, at the 
present stage, it is difficult to speculate the scaling behaviour 
of the S-wave velocity on stress, due to the interplay between 
the evolution of the stress anisotropy and the inherent fab- 
ric anisotropy. Indeed, the latter parameter is not included 
in most of analytical models within the framework of the ef- 
fective medium theory. 



Correlation between induced elastic anisotropy and 
stress-field anisotropy 



depend on both the induced stress-field anisotropy and the 
inherent fabric of the granular sample. The protocols of the 
sample preparation used here, i.e. de-compaction and rain 
deposition, allow us to obtain two distinct granular fabrics 
which also differ in the packing density <p as shown in Fig. 
[2] and probably in the coordination number Z. As suggested 
by eq. (Q}, we may investigate the correlation between the in- 
duced stress anisotropy and the elastic anisotropy via the ratio 
of elastic components dropping thus the parameters <p and Z. 

Fig. [8^ displays the ratio of the compressional moduli 
C33/C11 versus the induced stress anisotropy cr zz /cr xx for two 
granular samples with different fabric (loose and dense). The 
experimental data agree well with the effective medium model 
(eq. Q]), demonstrating that the anisotropy of compressional 
moduli is principally determined by the stress field anisotropy 
and is much less sensitive to the fabric of the granular pack. 
Also, we depict the evolution of the shear modulus ratio 
C44/C66 as a function of cr zz /o- xx in Fig. [HJ?. These results 
confirm again the correlation of the elastic anisotropy with 
the stress-field anisotropy, though the anisotropy is less pro- 
nounced for the shear modulus than for the compressional 
one. In contrast to the latter, the shear moduli are however 
very sensitive to the granular fabric. For the loose pack with 
isotropic fabric, the ratio of the shear moduli is in good agree- 
ment with the effective medium theory model based on the 
assumption of the isotropic distribution of contact angle. 
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In the above sections, we have observed that the elastic 
moduli and consequently sound velocities (P- and S-waves) 



FIG. 8: Correlation between the induced elastic anisotropy and 
stress-field anisotropy in the granular packs of different fabric. 



8 



V. CONCLUSION 

In summary, we have studied the anisotropic elasticity of 
dry glass bead packs using ultrasonic measurements. Both the 
influences of the fabric anisotropy and the stress anisotropy 
are investigated thanks to the granular samples of different 
inherent fabric under uniaxial loading in an cedometric cell. 
The stress anisotropy is characterized from the horizontal 
and vertical force measurements, whereas the elastic mod- 
uli are determined from velocity measurements of P- and S- 
waves propagating along the vertical and horizontal direc- 
tions, respectively. For the range of the applied stress in this 
work, our results show clearly that the compressional modu- 
lus anisotropy is more sensitive to the stress field anisotropy 
while the shear modulus anisotropy is more sensitive to the 
fabric anisotropy. 

We have tested the applicability of the effective medium 
theory to sound propagation associated with the incremental 
deformation. The predictions by the analytical model based 
on the affine approximation agree well with the measurements 



of sound velocity. Moreover, our velocity measurements of 
compressional waves confirm the scaling behaviour on stress, 
Vp ~ P 1 ^, where P is the normal component of the stress 
along the propagation direction. As for the shear waves, the 
stress dependences of velocities are more complex, being also 
strongly influenced by the fabric anisotropy which is not con- 
sidered in existing analytical EMT models. 

The present effective medium model is also applied to an- 
alyze the finite deformation of the granular pack. Our stress- 
strain measurements show the breakdown of this model to 
describe adequately the deformed state at large deformation. 
However, such discrepancy between theory and experiment 
can significantly be reduced when the friction between grains 
is switched off. As indicated previously lR l33ll . our observa- 
tion confirms that the model based on the affine approximation 
fails because it prevents the relaxation of the tangential stress 
at the grain contact, especially for large deformation. Other 
effective medium models are required for describing the large 
deformations of granular packs la, ll2[] . 
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